Effects of the magnetic field on the spaUation reaction implemented by BUU coupled 

with a phase-space coalescence afterburner 
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Based on the Boltzmann-Uehling-Uhlenbeck (BUU) transport model coupled with a phase-space 
coalescence afterburner, spallation reaction of Au at the incident beam energy of E,beam= 800 

MeV/nucleon is studied. We find that the number of test particles per nucleon has minor effects on 
the neutron to proton ratio (n/p) of the produced heavier fragments while it affects much on their 
yields. The external strong magnetic field affects the production of heavier fragments much than 
the n/p of produced fragments. The n/p of free nucleons is greatly affected by the strong magnetic 
field, especially for the nucleons with lower energies. 
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There has been a renewed interest in the study of 
spallation reactions induced by either nucleons or light 
charged nuclei, not only nuclear physicists but also as- 
trophysicists and nuclear engineers [ll, 0] . The spallation 
reaction is a kind of nuclear reaction in which a particle 
(e.g. proton) interacts with a target nucleus. Giving a 
high energy to the incident proton, the nucleus is then in 
an excited state and can de-excite by evaporation and/or 
fission. And then the high number of secondary neutrons 
are produced. The fragmentation can be used as a source 
of neutrons, for example, the ADS (Accelerator-Driven 
System) , which is as external source to drive the sub- 
critical reactor. The optimum proton energy for the pro- 
duction of neutrons by spallation in a heavy metal target, 
in terms of costs, target heating, and system efficiency, 
lies in the range 600 to 1000 MeV [1]. The spallation 
reaction is also related to the Spallation Neutron Source 
(a new generation of pulsed neutron sources based on 
the proton accelerator, and can provide high flux pulsed 
neutron beams for neutron scattering.), which is used to 
probe the structure of the microscopic world 5]. 

Proton induced reactions can be described by the semi- 
classical Boltzmann-Uehling-Uhlenbeck (BUU) transport 
model @, Q, which is quite successful in characteriz- 
ing dynamical evolution of nuclear collision. However, 
the description of fragment formation is not a trivial 
task. Because the BUU transport models are incapable 
of forming dynamically realistic nuclear fragments, cer- 
tain types of afterburners, such as statistical and coales- 
cence models, can be used as complements. The hybrid 
model can be usually used to study fragment formation, 
for instance, BUU-I- statistical multi-fragmentation mode 
(SMM) [^, BUU-I- nuclear multi- fragmentation model 
[9l-[l3|. However, these models all have their own short- 
comings, such as the lack of describing the dynamical pro- 
cess of the fragment formation or not giving the accurate 
excitation energies of fragments (which is very important 



for the de-excitation of hot fragment [T3| ) . 

In this work, we study fragment formation in proton- 
induced reaction p -i-^^'^ Au at the incident beam en- 
ergy of Eteam= 800 McV/nucleou using the BUU trans- 
port model coupled with a phase-space coalescence after- 
burner [13] . We first deduce a reasonable number of test 
particles per nucleon by fitting the experimental data. 
Then we discuss the phenomenon of multi-fragmentation 
in the external magnetic field and show how the magnetic 
field affects the spallation reaction. 

It is well known that the semiclassical Boltzmann- 



Uehling-Uhlenbeck (BUU) transport model [15| is quite 
successful in describing the evolution of nuclear collision. 
The BUU equation describes time evolution of the single 
particle phase space distribution function f{r,p,t), the 
equation is 



dl 

dt 



-:t+v- Vrf - VrU ■ Vpf = 7, 



collision 1 



(1) 



where f{r,p,t) can be viewed semi-classically as the 
probability of finding a single particle at time t with the 
momentum p at the position r. The right-hand side in- 
cludes the Pauli blocking. Because one has to initialize 
nucleon positions to start the cascade model, nucleon po- 
sitions imply a certain density. Then we can define the 
mean field potential U as a function of density. In our 
work, we use the Skyrme parametrization [Tsj . it reads 



U{p)^A{^) + B{^y 
Po Po 



(2) 
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Where cr — 7/6, A — -0.356 GeV is attractive and B = 
0.303 GeV is repulsive. With these choices, the ground- 
state compressibility coefficient of nuclear matter K= 201 
MeV. p = pn + Pp is the baryon density and p„, pp are 
the neutron and proton densities, respectively. 

The variations of nucleonic momentum are generally 
due to momentum and spatial location dependence of 
its mean-field-potential U, decided by the gradient force 
S/rU in Eq. (1). Besides the gradient and Coulomb forces 
added on the charged particles, the Lorentz force can 
also change the momentum of charged particle. For the 
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external magnetic field force of charged nucleon felt, we 
employ the Lorentz force equation [16|. According to 
Hamilton's equations, the propagations of nucleon are 

dpi/dt — —VrU{ri)+Fcoulomb + FLorentz, 

dri/dt = pj^m? -^pf. (3) 

The test particle method was first applied to the nu- 
clear Vlasov equation by Wong In the practical 
BUU calculations, to analyze fragment formation in pro- 
ton induced spallation reactions, we have used different 
numbers of test particles per nucleon. From them we 
choose 5 and 10 test particles per nucleon as reasonable 
values by comparing with experimental data. And in the 
next calculations to analyze fragmentation in the reac- 
tions with strong magnetic field, we use 8 test particles 
per nucleon as the reasonable value. The number of test 
particles per nucleon can not be too small as statistical 
fluctuations, which is inherent in any Monte-Carlo cal- 
culation, become apparent. Also the number can not 
be too large since we need to sample the phase space 
distribution of nucleons [11]. Usually we use 200~ 500 
test particles per nucleon to do BUU simulations, such as 
nucleonic emission or meson production. Different num- 
bers of test particles per nucleon really less affect BUU 
results. However, in the studies here, number of test par- 
ticles per nucleon affects much on the heavier fragment 
production. This is because small number of test par- 
ticles per nucleon gets more statistical fluctuations back 
when implemented by BUU transport model. While large 
number of test particles per nucleon smooths out statis- 
tical fluctuations (which cause fragment production). In 
fact, the number of produced heavier fragments is oftern 
very small compared with free nucleons. 

Because most BUU-type transport models are not able 
to form dynamically realistic fragments, some types of 
afterburners, such as statistical and coalescence models, 
are naturally used as a supplement to describe nuclear 
multi-fragmentation. Such kind of hybrid models can be 
used well to study fragment formation, for example, in 
the studies of collective flows of light fragments [Tj, \l§t - 
120( 1 and fragment formation in proton induced reactions 
Q. There are also some interesting works available on 
exerting advanced coalescence models [2Tl - [2^ . But these 
methods are in difficulties on predicting heavier frag- 
ments. For the purpose of our investigation, we take 
the simplest phase-space coalescence model [ij, [l^, [13] , 
a physical fragment is formed mainly when nucleons with 
relative momenta smaller than Pq and relative distances 
smaller than i?o- The results presented in the follow- 
ing discussions are obtained with Pq ~ 263 MeV/c and 
i?o = 3 fm. We have changed the two parameters prop- 
erly and find that they do not affect our studied here 
evidently. 

For the spallation reactions, the mass distribution of 
fragments are important for practical applications. Fig.[T] 
shows mass distribution of fragments in the p-\-^'^'' Au re- 
action at the incident beam energy of 800 MeV/nucleon, 



10' 
10^ 



c 

I 10^ 



10 



' I ' r 
p+Au@800 MeV 




10 test particles 



50 100 150 
Mass Number A 



200 



FIG. 1: (Color online) Mass distribution of fragments in the 
p Au reaction at the incident energy of ^beam~ 800 
MeV/nucleon. Theoretical calculations (with, respectively, 5 
and 10 test particles per nucleon) are compared with existing 
experimental data. Data are taken from [3]. 



the numbers of test particles per nucleon used in the cal- 
culations are 5 and 10, respectively. It is clearly shown 
from Fig. [T] that both results are roughly in agreements 
with the shape of the experimental data, especially in the 
zone of large mass number fragments {A > 130). In the 
mediate range (40 < A < 130), our results are not well 
in agreements with the experimental data. The results 
with 5 test particles per nucleon overestimate the produc- 
tion of fragments while the results with 10 test particles 
per nucleon underestimate the production of fragments. 
From physical point of view, small number of test par- 
ticles per nucleon gets more statistical fluctuations back 
when implemented by BUU transport model [l5| . Thus 
more fragments are produced. 

As discussed above, different numbers of test parti- 
cles per nucleon really less affect the BUU results. In 
the studies here, however, number of test particles per 
nucleon affects much on the fragment production (espe- 
cially for 20 < A < 150). The small number of test 
particles per nucleon gets more statistical fluctuations 
back when implemented by BUU transport model [15| 
while large number of test particles per nucleon smooth 
out statistical fluctuations. In fact, number of produced 
fragments of 20 < A < 150 is very small compared with 
free nucleons. The present studies are the first attempts 
to describe the produced heavier fragments. Further im- 
provements of our method are in progress. Because the 
experimental data are almost in the range of our simu- 
lations, we in the next studies take 8 test particles per 
nucleon as our model parameter. 

Fig. [5] shows neutron to proton ratio of nuclear frag- 
ments with, respectively, 5 and 10 test particles per nu- 
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FIG. 2: (Color online) Neutron to proton ratio of nuclear frag- 
ments with, respectively, 5 and 10 test particles per nucleon. 



cleon. We can see that the number of test particles per 
nucleon has minor effect on the neutron to proton ratio 
of nuclear fragments. The effects of the number of test 
particles per nucleon (with 5 and 10 test particles per 
nucleon, respectively) on the neutron to proton ratio of 
nuclear fragments are not more than 5%. 
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FIG. 3: (Color online) Mass distribution of fragments in the 
p -1-1^^ Au reaction at the incident energy of Ei,eam= 800 
MeV/nucleon with and without magnetic field, the strength 
of magnetic field added is 10^^ Tesla. 

The condition of strong magnetic field may exist in the 
universe, such as white dwarfs, neutron stars, and accre- 
tion disks around black holes, and the maximum value of 
magnetic fields in the universe may reach 10^° ^ 10"'^ G 
25 1 . And with the rapid development of laser technology 
obtaining strong magnetic field greater than 10 



Tesla artificially in terrestrial laboratory is possible. Also 
the strong magnetic field greater than SxlO^'' Tesla can 
be p rovided via energetic heavy-ion collisions technically 
|28l - l3l| . In the following studies, we set the strength of 
external magnetic field added to be 10^^ Tesla with the 
direction perpendicular to the reaction plane. 

To see the effects of external strong magnetic field on 
multi-fragmentation, we plot Fig. [3l mass distribution of 
fragments in the p+^^'^ Au reaction at the incident energy 
of Eiieam= 800 MeV/nuclcon with and without magnetic 
field. We use 8 test particles per nucleon in our simula- 
tions. From Fig. [Sj we can see that the magnetic field 
decreases most of the formations of fragment. Fragments 
are formed through statistical fluctuation in nuclear col- 
lision, the external magnetic field prevents protons to 
form cluster with other nucleons due to the Lorentz force 
added. However, as shown in Fig. 21 we can see that 
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FIG. 4: (Color online) Neutron to proton ratio of nuclear 
fragments in the p +^^'^ Au reaction at the incident energy of 
Ebeam^ 800 MeV/uucleon with and without strong magnetic 
field. 

effects of the external magnetic field on the n/p of pro- 
duced fragments are not very evident. 

The energy spectrum of free neutron is important for 
ADS and Spallation Neutron Source. We thus studied 
the effects of magnetic field on the free neutron produc- 
tion. Fig. [5] shows the energy spectrum of free neutron 
in nuclear spallation reaction p Au at the incident 
energy of Efceam= 800 MeV/nucleon with and without 
magnetic field. From Fig. [S] we can see that cross section 
of the free (especially energetic) neutrons produced in 
the spallation reaction with magnetic field is smaller than 
that without magnetic field. From Fig. [51 we can see that 
the neutron to proton ratio of free nucleons with mag- 
netic field is evidently larger than that without magnetic 
field. Free nucleons of the spallation reaction p -i-^^'' Au 
at ^beam= 800 MeV/uuclcon are mainly from the initial 
de-excitation of hot matter (executed by BUU through 
nuclear evolution) . Because the magnetic field holds pro- 
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FIG. 5: (Color online) Energy spectrum of free neutron in 
nuclear spallation reaction p+^^'' Au at the incident energy of 
Ebeam^ 800 MeV/uucleon with and without magnetic field. 
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FIG. 6: (Color online) Neutron to proton ratio of free nucleons 
in the p Au reaction at the incident energy of 'Ebeam — 
800 MeV/nucleon with and without strong magnetic field. 



tons to shoot out, the initial de-excitation of hot matter is 
partly prevented. We thus see small number of free neu- 



trons and protons in the spallation reaction p+^^'^ Au and 
large n/p of emitted free nucleons with external magnetic 
field. 

Our present coalescence method does not include frag- 
ment de-excitation processes, the fragments are all pre- 
fragments (their binding energies may not equal to their 
ground state energies). These pre-fragments may have 
de-excitation although we construct these fragments at 
the "final stage" of nuclear reaction. As shown in Ref. 
[s^ l. fragment de-excitations really affect the n/p of nu- 
clear fragments. The effects of fragment de-excitations on 
n/p ratios are evidently larger than the effects of mag- 
netic field on fragments' n/p ratios (shown as Fig. 0] in 
the text). However, we can see that the effects of mag- 
netic field on free nucleons' n/p (shown as Fig. [5] in the 
text) are much larger than the effects of fragment de- 
excitations on n/p ratios. Since our studies here are just 
the effects of the magnetic field on the fragment produc- 
tion and the effects of external magnetic field on fragment 
de-excitation processes are also less unknown, we thus in 
our studies neglect fragment de-excitation processes. As 
discussed above, the de-excitation processes of nuclear 
fragments do not affect our physical results here. 

As shown in Ref. [3d|, in heavy- ion collisions there 
is practically no electromagnetic field effect on observ- 
ables. This is due to the mutual dynamical compensa- 
tion of transverse components as clearly demonstrated in 
Ref. 3l| . Such compensation should be absent in a pure 



magnetic field and negligible in proton-nucleus collisions. 
Thus the n/p ratio is expected to be unchanged in heavy- 
ion collisions but its change for proton-induced collisions 
in external magnetic field should be emerged. 

In conclusion, based on the transport model BUU cou- 
pled with a phase-space coalescence afterburner, p + Au 
spallation reaction at the incident energy of Ebeam = 800 
MeV/nucleon are studied. The external strong magnetic 
field added may affect the formation of fragments and the 
n/p of emitted free nucleons. The effects of the Lorentz 
force on the n/p of heavier fragments are less evident. 
These studies may be useful for nuclear physics, astro- 
physics and nuclear engineering. 
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